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Abstract: We investigate a lattice model for Euclidean quantum gravity based on discretiza- 
tion of the Palatini formulation of General Relativity. Using Monte Carlo simulation we show 
that while a naive approach fails to lead to a vacuum state consistent with the emergence of 
classical spacetime, this problem may be evaded if the lattice action is supplemented by an 
appropriate counter term. In this new model we find regions of the parameter space which 
admit a ground state which can be interpreted as (Euclidean) de Sitter space. 
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1. Introduction 

The problem of constructing a theory of gravity which is consistent with quantum mechanics 
has a long history. The non-renormalizability of quantum general relativity can be traced 
to the fact that the Newton constant carries the dimensions of length squared corresponding 
to an action linear in the curvature. Actions containing higher powers of the curvature are 
usually thought to lead to violations of unitarity and can at best be thought of only as effective 
field theories |||]. 

Alternative approaches to gravity have tried to highlight the similarities to gauge theory, 
where in the case of gravity, the local symmetry corresponds to Lorentz invariance. In this case 
the corresponding gauge field is called a spin connection. It is necessary to introduce such an 
object when discussing fermions in curved spacetime where it is partnered by a new field the 
so-called vierbein which describes the local Lorentz frame and from which the usual metric 
tensor can be reconstructed. An action coupling the vierbein to the Yang-Mills curvature 
associated with this spin connection can be written down and shown to reproduce the usual 
Einstein equations under certain conditions - that the vierbein, regarded as a matrix, be 
invertible and that the torsion - the antisymmetrized covariant derivative of the vierbein, be 
zero. It is a first order formulation as both the vierbein and spin connection are to be varied 
independently to determine the classical equations of motion. 

This recasting of gravity in the language of gauge theory is termed the Palatini or Palatini- 
tetrad formulation of General Relativity Q . In principle this approach gives a natural starting 
point for a non-perturbative study of quantum gravity since it is possible to discretize the 
theory while maintaining exact gauge invariance using techniques similar to those employed 
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for (lattice) QCD. As such it is complementary to lattice approaches based on Euclidean 
dynamical triangulations see eg. |3|, ^ |5|, |6| and causal dynamical triangulations [0, ^, ^ and 
references therein. 

A number of proposals have been made earlier along these lines fl^, 0, |l2|, 13, 14] 



and numerical studies performed |15, 17]. Unfortunately this previous work encountered 
several difficulties; chief among these was the observation that the expectation value of the 
vierbein vanished in the classical weak coupling limit. This condition by itself is sufficient to 
invalidate the connection between the gauge theory and gravity as we will discuss in the next 
section. 

In this work we have explored the phase structure of a lattice model in which the discrete 
Palatini action is supplemented by an additional Yang-Mills term. We find a region of the 
parameter space of the model where the vierbein is non-zero, the torsion small and the 
curvature constant. 

The paper starts with a review of the Palatini formalism in the continuum, and describes 
its discretization on the lattice. The non-perturbative structure of this theory is then explored 
using Monte Carlo simulation and we show the classical ground state of the theory corresponds 
to vanishing vierbein and a curvature which approaches its kinematic limits ^. Thus the 
resulting lattice theory does not possess an appropriate ground state which can be identified 
with a continuum geometry. We then modify the model to incorporate an additional curvature 
squared operator and survey the expanded phase diagram finding indications of a region 
where the ground state of theory can be identified with Euclidean de Sitter space. In the 
final section of the paper we discuss what these results imply and what further work must be 
done to solidify any possible connection to a theory of gravity. 



2. Review of the Palatini formalism 

The Palatini action can be regarded as the most general action constructed from the curvature 
of the spin connection and the vierbein which is invariant under local Lorentz transformations 
and whose definition is independent of the choice of background metric. It takes the form 

^Palatini = ^ Jd'x e,,xp e'^^' e\ ei {r^^ -^e\ e^,) (2.1) 

with Ip = \/AttG the Planck length. The usual metric of General Relativity is then related to 
the vierbein via 

9,iu = r]ije'^ei (2.2) 

where rjij = 6ij is a flat (Euclidean) metric which henceforth we shall leave implicit in formu- 
lae. This relation is clearly invariant under local 5*0(4) gauge transformations mediated by 

^In Euclidean space we replace the Lorentz group SO{3, 1) by its compact analog 50(4). In the context of 
the lattice theory this ensures that the local curvature is bounded both above and below 
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the spin connection cj^ where 



u, = Y,up'^ i<j = l...A (2.3) 



with T^^ being appropriate generators of 50(4). The curvature R]iu{^) appearing in eqn. 2.1 
is the usual Yang-Mills field strength associated with the spin connection. 

Variation of this action with respect to the spin connection yields the torsion free constraint 
^[fj. ^V] ~ ^ which can be used to solve for the spin connection in terms of the vierbein provided 
the inverse exists. The solution is 

where, 

%up = {dpet-d,e';,)epk (2.6) 
Variation of the action with respect to the vierbein then yields the equation, 



e'^'''e,,,,{E^,--^ei)e>l = (2.7) 



A 
3" 

which, in conjunction with the torsion free solution given above uj = a;^(e), just reproduces the 
usual Einstein equation for pure gravity with a cosmological constant. The vacuum solutions 
are then constant curvature spaces - here Euclidean de Sitter space - the sphere S^. 

As MacDowell and Mansouri noted |]l^] it is possible to forge an even stronger connection 
to Yang-Mills theory by combining the 5*0(4) spin connection and vierbein into a single 
gauge field associated with a enlarged 50(5) gauge symmetry 



= u'^T'^ + -e'^T''' i,j = l,..A (21 



1 

— I 

/ 

where / is some scale inserted to render the vierbein dimensionless. In a similar fashion 
the 50(5) curvature F^^, can be decomposed into components transforming under the same 
50(4) subgroup 

Fpv — [Rp,u ~ 'j2^[fJ-^v\ j ^ 

+ yl)[^etjr5' i,i = 1...4 (2.9) 

with R the 50(4) curvature. 

An 50(4) invariant action quadratic in this 50(5) curvature may then be written down 

S = Kj d^xe,,,pe^^'''Fi^,F^j, (2.10) 
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By expanding the 50(5) field strengths according to eqn. |2.9| this is easily seen to be nothing 

more than the Palatini action given in eqn. (|2.1j ) with the dimensionless coupling k = ^ (^j-^ 

giving the size of the Universe in units of the Planck length and A = Notice that the 
classical equations of motion derived from the Palatini action just correspond to setting the 
SO (5) curvature to zero. 

This action while formulated in terms of an 50(5) connection exhibits an explicit breaking 
to 50(4). It is natural to generalize this action slightly to try and realize this as a spontaneous 
breaking. The following 50(5) action accomplishes this at the expense of introducing an 
additional scalar field (j). 

S = I d'xe.^^^pe'^'^'^FliF^lr (2.11) 

On the assumption that the scalar acquires a vacuum expectation value of the form (f)^ = k6^"^ 
we recover the previous action. 

The appearance of a scalar field (j) ^-nd an associated 50(5) gauge symmetry seems at 
first sight to be somewhat arbitrary. However, it is possible to show that this 50(5) invariant 
theory arises naturally as a compactification of a Chern-Simons theory in five dimensions. 
This Chern-Simons theory is a generalization to five dimensions of Witten's formulation of 
three dimensional gravity as a Chern Simons theory [?, ?]. In the Euclidean case considered 
here this amounts to constructing a topological gravity theory in five dimensions with internal 
symmetry group 50(6) and action 

/" 3 3 

5cs=/ <^7A7^A7^ + -^]A^]A^lA7^ + -^lA^lA^2A^7A^l> (2.12) 
Jms 2 5 

where for brevity we used the language of differential forms to represent the contraction of 
spacetime indices with the five dimensional epsilon symbol. The Lie algebra valued connection 
is 

%=J2^^^Jab;A,B = 1,...,6; (2.13) 

A<B 

with Jab = [7A)7s] the generators of 50(6) and TZ the curvature. The the angular brackets 
indicate that the group indices are contracted with the 50(6) invariant tensor ^abcdef- This 
Chern-Simons Lagrangian is related to the Euler density in six dimensions via the relation, 

dies = eABCDEF n^'^ A 71^^ A (2.14) 

Stokes' theorem then guarantees that the Chern-Simon's theory is invariant under local gauge 
transformations up to possible boundary terms. 

It is useful at this point to decompose the connection in terms of quantities transforming 
simply under an 50(5) subgroup 

J] = ^A'^^J,fe + ^SV6a a, 6 = 1... 5 (2.15) 

a<h a 

term quadratic in the 5*0(4) curvatures also appears / e^ji/Ape'"'''' J?Jf^i?Ap but can be neglected as it 
corresponds to a topological invariant - the Euler number 
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Similarly the SO (6) field strength decomposes according to 



n = J2(F"'' + E^^E'')Jab + J2DE^Jea (2.16) 

a<b a 

where F is the SO (5) curvature and DE the torsion. The action up to boundary terms 
becomes 

Scs = f eabcde ( A F^^ AF'^'' + -E'' AE'' AE^ A F'^^ + - E" A E^ A E^ A E'^ A E' 
Jms \ 3 5 

(2.17) 

In this form it is clear that the action does indeed describe a gravitational theory contain- 
ing a "Gauss-Bonnet form" , an Einstein term and a cosmological constant in five dimensions. 
Furthermore, this theory is invariant under the full (Euclidean) de Sitter group not just local 
Lorentz transformations. 

To make contact with the MacDowell-Mansouri form of the Palatini action we must 
compactify this theory down to four dimensions. Let us assume that the five dimensional 
manifold is of the form M5 = M4 x S^ /Z,2 with Z5 the extent of the fifth dimension. On 
the four dimensional boundaries we will assume that the gauge field satisfies the following 
boundary conditions, 

= 76 ri^76 ; /i = 1, . . . ,4 (2.18a) 
^5 = -16^516 (2.18b) 

This breaks the 5*0(6) gauge symmetry down to 50(5) at the boundaries with the only 
surviving four dimensional fields being 

^5, <• (2.19) 



Returning to the action in eqn. 2.17 it shoul be clear that only the "Gauss-Bonnet" term 



survives at tree level on the four dimensional boundaries. Indeed, this term becomes 

J d^x eabcde e^^Ap El F'^l Ft; (2.20) 

This is nothing more than the previous action with the fifth component of the five dimensional 
vielbein playing the role of the scalar field in the four dimensional theory. The development 
of a vacuum expectation value for the scalar field could then be realized in terms of the 
appearance of a non vanishing Polyakov line in the fifth direction extending between the two 
four dimensional boundaries. 



Thus the S0(5) invariant theory given in eqn. 2.11 can be realized by an appropriate 
compactification of a topological gravity theory in five dimensions. The possible existence 
of an underlying exact 50(5) gauge symmetry in the four dimensional theory has some 
advantages; it allows us to restrict possible counter terms to those invariant under 50(5), it 
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makes the choice of the Haar measure for 50(5) natural in the path integral and also, as we 
now explain, it makes the restoration of diffeomorphism invariance a more likely possibility. 

In [?] Witten showed that in fact diffeomorphism invariance is intimately connected to 
invariance under both lorentz transformations and local translations. In the four dimensional 
model considered here the result of a general coordinate transformation with parameter —^'^ 
acting on the spin connection and vierbein can be written compactly as a gauge transformation 
with parameter £,^A^ acting on the 50(5) gauge field A = u -\- e plus a term which vanishes 
on flat connections with F = 

5^A^ = -D^ii'^A,) - ^F^, (2.21) 

Thus provided the theory is S0{5) invariant and we consider only small fluctuations around 
a flat background the theory will automatically be invariant under general coordinate trans- 
formations. 



3. Lattice theory 



It is straightforward to formulate the Palatini action in the form given by eqn. 2.10 on a 
hypercubic lattice as was first suggested in ||Tot| . 

S'p = K ^ ^ e^^ApTr {-if>U^uUxp) (3.1) 

where 

U^. = U^{x)U,{x + fi)Ul{x + u)Ul{x) (3.2) 

is a Wilson plaquette variable and takes its values in the group SO{5f. Notice that m our 
study we have taken as generators of SO{5) the matrices T*-' = ^[7*, 7-'] where 7^, z = 1 . . . 5 
are the usual four dimensional Dirac matrices together with the chiral matrix 75. This differs 



from the earlier numerical work reported in [16, 17 1 which utilized the vector representation 



of SO (5) and implies that we are actually simulating a lattice theory based on the covering 
group spin(5). As we will see this will turn out to be rather important. 

Notice though that the introduction of a lattice has necessarily broken the coordinate 
invariance of the theory in the base space. Thus one should expect that quantum corrections 
will generate additional operators whose structure depends on the existence of this background 
lattice. Assuming that these operators depend only on the S0{5) curvature there are 3 such 
terms in addition to the Palatini action which are (marginally) relevant by power counting; 

• / d^xe^,,ApTr {F^^Fxp) 

. fd^xTi{F^,F^,) 



^For simplicity in this study we have not symmetrized our action to ensure rotational invariance nor have 
we concerned ourselves at this point with issues of reflection positivity 
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fit/ j 



The first of these is a topological invariant in the continuum analogous to the instanton 
number and hence can be neglected. In principle we should include lattice versions of both 
the second and third terms. In practice we have focused these initial investigations on just 
the second term. On the lattice we implement it with the usual Wilson plaquette action 



(3.3) 



Our choice of this term has clear motivation; it clearly will favor configurations in the large 
K limit with vanishing 5*0 (5) curvature - a feature which will see shortly is not true of the 
pure Palatini action. Such configurations are minimally required to achieve a connection to 
classical gravity. It also has the merit of removing potential lattice doubler modes which will 
be present in the Palatini action as was first observed in p^ . 

Of course the presence of a Wilson term is not compatible with coordinate invariance in 
the base space. It will clearly be very important to test for a restoration of this property in any 
continuum limit. Nevertheless the important point to realize is that this term is necessarily 
induced in the lattice theory and to have a hope to obtaining the correct continuum limit we 
should include it in the bare lattice action and tune its coupling appropriately as the lattice 
spacing is reduced. The lattice action we consider then has the form 



S = K {aSw + (1 — Oi)Sp) 



(3.4) 



which allows us to interpolate between the pure Palatini action defined by a = and pure 
Wilson action when a = 1. 

Finally, we need to give a presecription for extracting the various 5*0(4) components of 
the connection and curvature from the basic SO (5) variables in the theory. We used the 
simple expressions for the vierbein and torsion 



= Tr [T^^C/^ 



^ fit/ 



Tr 



^ 2 (^A**^ ^i^fJ-] 



while the SO{4) components of the curvature TZ^u = R^iu — are given by 



(3.5) 
(3.6) 



= Tr 



i = 1...4 



(3.7) 



As in the previous studies [O, 16 1 we have assumed an 5*0(5) invariant Haar measure on 



the group in the path integral defining the quantum theory. We have employed a standard 
metropolis algorithm to perform the Monte Carlo simulation. 
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Figure 1: Expectation value of Palatini action for a = vs coupling k 

4. Numerical results 
4.1 Pure Palatini 

In fig. |l] we show the expectation value of the pure Palatini action (a = 0) as a function of 
the coupling k for a sequence of three lattice sizes L = 6, 8, 12. One sees that the action 
approaches its lower kinematic bound of Sp = —24 (which refiects the number of non-zeroes 
of the e symbol) independent of L as the weak coupling limit k — )■ oo is taken. This limit is 
reached when the plaquette variables satisfy the extremality condition 

U^,u = -Ib^^^vpxUl^ (4.1) 

We will see later that pure Palatini configurations satisfying this condition possess a maximal 
50(4) curvature and hence cannot be interpreted as corresponding to a smooth spacetime. 
Presumably this refiects the usual unboundedness problem of Euclidean quantum gravity. 
Furthermore, we observe that the expectation value of the vierbein approaches zero as k — >• oo 
as can be seen in fig. ^ in which Ti{g) = ^ e^e^ is plotted as a function of k for the same 
range of lattice sizes^. As emphasized earlier the vanishing of the vierbein implies that the 
corresponding metric is zero and removes the possibility of interpreting the classical 
equations of motion as corresponding to the field equations of a metric theory of gravity. 

Before discussing the situation with a > we first return to an issue of what represen- 
tation should be used for implementing the local Lorentz invariance. As we have described 

*We absorb the scale / into our lattice vierbein throughout this paper 
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in the last section, all the results shown in this paper have been obtained using the four di- 
mensional representation of spin(5) given in terms of commutators of four dimensional Dirac 
matrices. As one can see in the previous two plots this choice generates a smooth dependence 
of observables on the coupling k. This is quite different to what was seen in the only previous 
numerical study of this system where a very strong first order phase transition was observed 
separating strong from weak coupling.^ Indeed the observed hysteresis effects were so severe 
in the latter case that it was hard to even thermalize the lattices in the weak coupling phase. 
This earlier study utilized the fundamental or vector representation of SO{f)). This observed 
disparity in the phase structure of these two lattice theories depending on choice of represen- 
tation is seen even in the pure Wilson theory (a = 1) as can be seen in fig. |3| which shows 
the average action plotted as a function of n for the two representations While the two repre- 
sentations agree at weak coupling as they must since they possess the same Lie algebra, they 
differ at strong coupling and the lattice action employing the fundamental representation of 
SO{b) suffers a strong first order transition for k ~ 0.75 which is completely absent in the 
spin(5) representation. This is analogous to the situation for the lattice theories of 50(3) and 
SU (2) - the former exhibiting a first order bulk phase transition separating weak from strong 
coupling while the latter does not. Indeed, the analogy is even stronger since SU{2) = spin{3) 
which leads us to conjecture that while the SO{N) groups have first order bulk transitions 
their covering groups spin{N) do not. Certainly from a practical point of view the use of the 
spinor representation is clearly superior to that of the vector representation. 

^One needs to rescale k in this work by 16 to compare the two couplings in these simulations ||l^, FH| 
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Figure 3: Comparison of vector and spinor representations for 5*0(5) Wilson action 



4.2 Palatini plus Wilson 

In this section we shows results obtained in the model with a > 0. First consider the 
expectation value of the Palatini action in such models. Fig. ^ shows this for four different 
values of the coupling a = 0.0, 0.25, 0.5, 0.75 for a fixed 8^ lattice as a function of k. Clearly 
the curves fall into 2 classes; for small a the expectation value of the action approaches the 
kinematic boundary corresponding to maximal (negative) curvature while for larger a they 
approach vanishing action in the weak coupling limit. A similar contrasting behavior is seen 
in the plot of the vierbein as revealed in fig. |5|. At small a the vierbein is driven to zero for 
large k as for pure Palatini while for a above some critical value ot the vierbein remains 
non-zero in the weak coupling limit. Thus it appears that a sufficiently large coupling to the 
Wilson operator stabilizes the vierbein. However this is not enough - we also require a zero 
torsion condition and a small curvature. We now turn to these other observables. Fig. ^ 
shows a plot of TZ^iiuT^^iiu and TJ^^TJ^^^^ for the same values of a as a function of the coupling k. 
For small a the expectation values are driven to values similar to the case for pure Palatini - 
the 50(4) curvature attaining near maximal (negative) values as k — )• oo^. However, again, 
for large values of a both the torsion T and curvature TZ approach zero at weak coupling. 
Indeed, notice that for a = 0.75 the torsion and curvature are degenerate to within small 
errors in this limit. 

These are quite encouraging results but one might worry that the Palatini action is 
playing no role at all for large a and the physics is being dominated by just the Wilson term 

^Notice that 7^ ^> i? for e = 
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Figure 4: Expectation value of Palatini action vs coupling k for various a 



- which would invalidate any connection to gravity. To show that this is not the case we have 
plotted a "self-dual" order parameter given by^ 

On the extremal configurations given by eqn. [4.1| this attains a value of minus one. This 
should thus be the case for pure Palatini at large k. In contrast it should vanish in the case 
of the pure Wilson action. In fig. |^ we show the value of this order parameter for the same 
range of a and k on the 8^ lattice. Notice that the curves for a = 0.5,0.75 indicate that 
order parameter appears to approach a constant for large k at a value which is intermediate 
between the pure Palatini and Wilson values. We regard this as a piece of evidence disfavoring 
a system governed solely by the Wilson term. 

Up to this point we have concentrated on showing results as a function of k for several 
discrete values of a. However it is also very instructive to show the data as function of a 
at fixed k. For example, the expectation value of the Palatini action is shown in fig. ^ For 
a < ax with ax = 0.4 — 0.5 it is clear that the Palatini action fiows to its kinematic lower 
limit as K — )• CO independent of the value of a within this range while for a > ot the trend is 
reversed and the action goes to zero for large enough k. This suggests two different ground 



'^This is analogous to the self-dual order parameter introduuced in ||l^ although in that case the duahty 
operation involves only the internal space 
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Figure 5: Expectation value of TT{g) vs coupling k for various a 



states are available to the system at least for large enough k. whose domain of attraction is 
determined by the coupling a. For a < ax the system at weak coupling is governed by the 
pure Palatini action while for a > ot the system flows to a ground state corresponding to 
vanishing 5*0(5) curvature. In the limit a — )• 1 this is just the usual Wilson vacuum. To 
contrast with this behavior notice that for small k coupling eg k = 0.5, the value of the action 
depends only very weakly on the associated coupling a and a system seems to exist in just 
one phase. Of course one expects that the ground state of the system for small k is dominated 
by lattice artifacts so this latter phase is ultimately uninteresting. 

This cross-over between two different ground states for large k is clearly seen by looking 
at the torsion and curvature as a function of q. A plot (fig. |9|) showing T and TZ for k = 20.0 
as a function of a confirms the importance of this threshold coupling ay As a increases the 
curvature drops towards to zero rapidly merging with the torsion for a > ax ~ 0.5. 

A similar picture is revealed by looking at the plot of the order parameter P as shown 
in fig. 1^. It appears that the order parameter equals minus one for all a < ay as k — )• oo 
but then rises towards zero as a varies in the range a > ot- Notice also that there is little 
K dependence in the curves for a > ax - the curves for different k rapidly approach fixed 
well-defined envelope in this region of parameter space. Again, let us emphasize that the 
small K data doesn't fit this picture - there is no evidence for a threshold value of a in the 
K = 0.5 curve. Thus the phase diagram seems to contain a single strong coupling phase but 
the possibility of two dramatically different phases at large k coupling. It is unclear from 
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Figure 6: Curvature and torsion components of the 5*0(5) field strength vs k 
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Figure 7: Order parameter P as a function of k for several a 



the data whether there is a discontinuity or not for a — >• ot- One hint favoring a true phase 



transition at ot can be obtained by looking at the local volume element e*-''^'e^vApe* eie^e' 
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Figure 8: Expectation value of Palatini action as function of a for several k 
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Figure 9: Curvature and torsion as a function of a for k = 20.0 



We can extract a gauge invariant measure for this on the lattice by extracting an expression 
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Figure 10: Order parameter as a function of a for several k, 



for a matrix valued representation of the vierbein at a site: 

e^{x)=^^U^{x)^^Ul{x)-I (4.3) 

which leads to an expression for the gauge invariant volume element 

y/g = e/.i.ApTr {e^euexep) (4.4) 

The absolute value of this quantity is shown in fig. 11 for k = 20.0 and L = 8 as a function 
of a. For small a the vierbein vanishes and so does the volume element while for a > ot 
it becomes constant. However, in the vicinity of ot we see a sharp spike corresponding to a 
large local volume element in lattice units. 



5. Discussion 



We have simulated a Euclidean model for quantum gravity based on discretization of the 
Palatini action including a cosmological constant term and supplemented by a Wilson term. 
Unlike a previous study of the pure Palatini action we have employed the spinor representation 
of the spin(5) covering group to define the lattice action. This is faster to implement and 
appears to avoid some of the strong lattice artifacts seen in earlier studies of this action. 

We show that the pure Palatini lattice theory is sick suffering from a local curvature 
which attains the maximal (negative) value consistent with the compact gauge symmetry, a 
vanishing vierbein and lattice doublers. 
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Figure 11: Absolute value of the volume element vs a ior k = 20.0 



We argue that the Wilson operator serves both to regulate all of these problems and 
furthermore will necessarily be induced at the quantum level if it is not present in the bare 
lattice action. In the expanded parameter space determined by the Palatini and Wilson 
couplings we find evidence for two new ground states for sufficiently weak coupling (large k) 
— one corresponding to pure Palatini and suffering from all the old problems and a second in 
which the vierbein is non-zero, and the torsion and curvature are small. These two vacua can 
be distinguished by a self-dual order parameter and can be accessed by tuning the coupling 
a which controls the mixing with the Wilson term. 

It appears that for a couplings close to some threshold value ot the new vacuum is not 
just that of the pure Wilson action but retains a memory of the Palatini term as revealed 
by the self-dual order parameter acquiring a non-trivial value there. It remains to be seen 
whether this threshold value ar can be thought of as a true critical value and if so whether 
this critical value corresponds to a continuous or discontinuous phase transition. The latter 
issue is of course a crucial issue to address in the context of obtaining a non-trivial continuum 
limit. 
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